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ABSTRACT 



A first order phase transition leading to deconfinement and chiral restoration 
is a likely possibility for QCD, at least in some region of the temperature- 
density plane. A signal for a unique transition is that the order parameters 
for such transitions (which can be understood in terms of symmetries only in 
limiting situations of very massive or massless quarks) are both discontinu- 
ous at the same critical temperature. We show that such a situation can be 
understood on a precise thermodynamical basis because of a general relation 
among discontinuities which holds for first order transitions. We derive the re- 
sult by a generalization of the Clausius-Clapeyron equation and also through 
the effective action formalism. We illustrate the theorem in an elementary 
example. 



1 Introduction and summary of conclusions 



There are two limiting situations in QCD where one expects the occurrence of phase 
transitions characterized by corresponding symmetry breaking. In the hmit of infinite 
quark mass, the deconfinement transition can be related to the breaking of a discrete 
group Z{3), which is a symmetry of the QCD lagrangian in that limit. The Polyakov 
loop acts as order parameter for the transition. In the opposite limit, of massless quarks, 
the chiral phase transition is related to restoration of chiral symmetry. The fermionic 
condensate acts as the relevant order parameter. 

In the realistic case of QCD with quarks of finite mass, neither the discrete Z{3) group, 
nor the global chiral symmetry are good symmetries of the lagrangian. One may however 
still consider to look at what happens to the two preceding order parameters, though they 
had been constructed for differing limiting situations. 

A computer simulation, with 4 light quark flavors with dynamical mass m/T ^ 0.1, 
has shown an abrupt transition for the two order parameters at practically the same 
temperature 0, 0. 

Another important issue is the order of the transition. Recent results, reviewed in ref. 
101 , point towards a first order phase transition in the case of two almost massless quarks 
and a massive one. It should also be mentioned that there appears to be a sort of critical 
value for the mass of the third quark, above which the transition becomes second order. 

The problem of the transitions has also been studied in ref. [^] assuming an effective 
lagrangian. The authors argue that the deconfined phase should be characterized by the 
vanishing of the gluonic condensate, {0\F^^F'^^\0)P{g)/ {2g)=0, and furthermore that the 
deconfinement transition should be first order. For this discussion it is convenient to 
introduce separate critical temperatures and for the chiral and the deconfinement 
transitions, as one would expect when ignoring couplings between different degrees of 
freedom. The authors of ref. [^] show that, if > T^, the first order deconfinement 
transition drives a discontinuity of the fermionic condensate. As a consequence the chiral 
transition appears of first order, whichever was its nature when ignoring other degrees 
of freedom, and at the same temperature of the deconfinement one. We also mention 
indications that the topological susceptibility drives the transition of other quantities, 
such as the quark condensate or the Polyakov loop. 

We shall prove here a general thermodynamical result. We shall show that a gener- 
alized version of the Clausius-Clapeyron equation relates discontinuities of the thermal 
averages of different observables at the critical point of a first order phase transition. In a 
more precise way, we consider a general hamiltonian, written in the form H = J2a ^aC>a, 
where are constants and Oa arbitrary operators. The critical point of the first order 
phase transition will depend, in general, on a subset of the parameters A^, to be called 
{Afc}. Then all the thermal averages of the operators Ok, related to the A^'s, are expected 
to become discontinuous at the critical point, both in their dependence from temperature 
and from density. 

The consequence for QCD is obvious, assuming that there is a first order transition and 
that the critical point depends on parameters, such as quark masses, the QCD coupling 
constant, etc. Then the quark condensate, the gluon condensate, etc., are expected to 
exhibit simultaneous discontinuities at such a critical point. Of course, this is not a 
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proof that such a scenario necessarily occurs in QCD, but it clearly shows that such an 
apparently curious thing as a unique phase transition, perceived through order parameters 
related to different, even incompatible, limiting situations of the theory, has a simple and 
precise thermodynamical basis. 

For the derivation we first present a generalization of the Clausius-Clapeyron equation 
(section 2). In section 3 we deal with the discontinuities. Wc give a derivation based on the 
functional effective action formalism in section 4. A most elementary example, developed 
in section 5, illustrates the simplicity of the result. 



2 Generalization of the Clausius-Clapeyron relation 

To derive the relation we are interested in, let us start by recalling some thermodynamical 
relations. Consider the partition function in the Grand Canonical Ensemble 

Z^Tr e-^(^-'^^) (2.1) 

and a Hamiltonian of the form 

TV 

H = Ho+Yl Aa 0« (2.2) 

where are real parameters conjugate to the hermitean operators Oa- The expectation 
value of the operator Oa is 

id.) ^ - ^ 



OK 



da e-^^^-'^^^ (2.3) 



In the same way, by varying with respect to the chemical potential /x, one has the 
mean particle number (or it could be any conserved charge, and we are not necessarily 
considering a fermionic theory) 



(N) = (2.4) 



We also recall that the entropy S is 

d (T In Z) 

Now, let us assume that the theory describes a system which can exist in two phases, 
each described by an equation of state Pi = Pi{T, n) (i = 1, 2), where pj is the pressure 
in the particular phase; and that there exists a set of points (//c, in the plane of 
chemical potential and temperature, by crossing which the system undergoes first order 
phase transitions. In any of these points the equations of state of the two phases coincide 

Pi(Tc({A,}), //e({A,}), {Afe}) = P2(7;({A J), /Xe({Afe}), {Afe}) (2.6) 
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({Afc} indicates the subset of parameters which affect the critical behaviour). 

If the pressure depends smoothly on {A^}, we can require the equality of the total deriva- 
tives of the two members in eq. (2.6) with respect to 



dpi 



dp2 



obtaining 
dpi 



dXk 



+ 



dpi 

or 



dJc 
c dXk 



+ 



dpi 
dfjL 



dXk c dXk 



duc _ dp2 
" dXi 



dXk 



+ 



dp2 
dT 



dT^ 
c dXk 



+ 



dp2 



duc 
" dXk 



(2.7) 



(2.^ 



Here the subscript means that the quantities are evaluated at (/Xc, T'c) • We recall that 

c 

according to the Ehrenfest classification, in a first order phase transition the pressure 
(or any another state function) is continuous whereas its first (partial) derivatives are 
discontinuous. Also, in the thermodynamic limit 



p V = T In Z 



(2.9) 



(here V is the three-dimensional volume and the pressure has been taken vanishing at 
T = /X = 0) Thus eq. (2.8) can be written in the form 



disc (ofe) = — — ^ disc s + -j^ disc n 
dXk dXk 



duc 



where 



(Ok) 



s 



n 



(N) 



(2.10) 



(2.11) 



V V V 

Let us now imagine two lines of first order critical points, one at A^ and another at Xk+dXk- 
If we consider eq. (2.10) at {fic,Tc) on the line at A^, we can unambiguously evaluate the 
discontinuities of the densities (6^), s, and n, whereas we have to specify on which point 
{/j,^, T^) of the line at A^ -|- dX^ we end in order to calculate the two derivatives dTc/dX^ 
and d/ic/dXk- In other words, although these derivatives are bound to satisfy eq. (2.10), 
we have the freedom to choose the direction along which to evaluate them. In particular 
we can consider the variation of (of fic) with respect to A^ at fixed fi (at fixed T), and 
in this way we can express disc (6^) only in terms of of disc s (of disc n), eliminating the 
other discontinuity. Thus, under these conditions, we have the further relation 



dT,. 



dX, 



disc s 



dfic 



dX, 



T 



disc n 



(2.12) 



The previous result can be also expressed by the fact that, in presence of a critical line, Tc 
and are not independent variables and that if we fix for instance the chemical potential 
dTjdXu = dTjdXu 

Let us now consider the thermodynamic relation 



rr, 9p dp 

e = -p + r J 

= — p + Ts + fin 
which allows to determine the latent heat for a first order phase transition 

disc e — Tf. disc s -\- disc n 



(2.13) 



(2.14) 
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By using eqs. (2.10) and (2.12), disc e can be entirely expressed in terms of disc (6^) 



/ 



disc e 



1 

d{\YiT,)ldh 



+ 



9(ln/ic)/9Aj 



disc (dfc) 



(2.15) 



The relation (2.15) is of the type of the Clausius-Clapeyron equation for a liquid- vapor 
first order transition, 

dTc disc V 



disc s 

[v is the specific volume) as shown in a recent paper by H.Leutwyler 
potential. 



for zero chemical 



However the relation (2.15) is quite more general; it concerns any operator and it 
allows for a chemical potential. We must remember that to get (2.15) in this precise 
form, eq. (2.12) is necessary, which requires the existence of lines of first order critical 
points at varying parameters {A}, rather than an isolated critical point. 



3 Relations among discontinuities 

Eqs. (2.10), (2.12) (and (2.15)) hold for the thermal average of any operator whose 
conjugate parameter influences the critical point. If Oi and Oj are two of such operators, 
eliminating alternatively disc s and disc n from eqs. (2.10) and (2.12), written by varying 
with respect to the parameters Aj and Xj, we also get 

^"'^ = H/fl: ^"^^ ^'-'^ 
^^^^ ^'^^ = w!m^T "^'^ ^'^^ ^'-'^ 

Namely, all the operators conjugate to parameters which affect the critical point are 
expected to be discontinuous at the phase transition point. The ratio of their discontinu- 
ities is simply the ratio of the variation of the critical temperature or chemical potential 
with respect to the correspondent conjugate parameters. From the two previous equations 
one also gets 

dTjdXi\^ ^ dfiJdXilT , . 

dTjdX.l ~ diiJdX,\T ^ ■ ^ 

Of course we could have employed other variables instead of the pressure and followed 
the procedure starting from eq. (2.6) obtaining relations similar to eq. (2.15). For 
instance, by equating the chemical potentials /ij(p, T) in the two phases for the liquid- 
vapor transition, varying with respect to p {Tc is a function of p), one obtains eq. (2.16). 
The reason why we have chosen the pressure will be clearer in the following. 

In this context our aim is to stress the essential ingredients. If in some system a first 
order phase transition occurs, using the fact that the state functions are continuous, but 
not their partial derivatives, one can easily relate the various discontinuities appearing in 
the description, obtaining alternative expressions for quantities such as the latent heat. 
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In particular these relations can be useful, as already remarked by Leutwyler [Q, if one 
can easier compute the discontinuity of some operator rather than the discontinuities of 
the derivatives of the pressure. 

The previous analysis can readily be extended to cases where there is a set of chemical 
potentials /x, coupled to charge operators iVj. 



4 The effective action derivation 

Let us consider a field theory at finite temperature T and chemical potential /i, described 
by an euclidean effective action F. The effective action is a functional of a set of operators 
{A}, the dynamical variables of the theory, 



r = r 



{A}-{a}-T,fi 



(4.1) 



In the previous equation we have also explicitly written the dependence of F on a set of 
parameters {a}, which are external sources conjugate to the operators {A}. Suppose that 
by varying the temperature T and/or the chemical potential /i, the effective action can 
describe different phases of the theory, and that there exist lines of first order critical points 
[Tc, fic) in the (T,fj,) plane. Then in any of these points (Tc,^c), T has two degenerate 
absolute minima {A}2, where we can write the condition 



r {i}i;{a};T,({a}),/i,({a}) =T {ijs; W; T,({a}), /!,({«}) 



(4.2) 



Then we equate the total derivatives with respect to any parameter ak which affects the 
critical behaviour, getting, apart from terms continuous at the transition. 



dT 



+ 



dT dTr 



+ 



dT 



+ 



dT d{A} 



dak dT dak dfidak d{A} dak 



Mil 



same 



{A}2 



(4.3) 



The derivatives of T with respect to {A} vanish on both sides of the previous expression, 
since they are evaluated at the two minima of the effective action. Actually the effective 
action evaluated at these points is nothing but — InZ and therefore — p V/T (since it 
is defined for infinite volume) for the two coexisting phases. Thus at the physical point 
it looses its dependence on the dynamical variables, whereas it remains a function of the 
sources {a}. Furthermore, in the case of constant fields, we can employ the euclidean 
effective potential 



V = = -p 

13V ^ 



(4.4) 



which directly gives the pressure. 



It follows that eq. (4.3), written for the effective potential, is equivalent to eq. (2.8). 
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5 A simple example 



Let us now consider a simple example. It is almost trivial, as it involves a single field. 
However it is useful to understand how the theorem works. 

Consider an effective potential containing powers up to the sixth order of a scalar field 
0, which plays the role of an order parameter [0: 

y = a(j)'^ - b(f)^ + c(j)^ + h(j) (5.1) 

where we assume c > and h small. We will show how the equations (2.10), (2.15) are 
satisfied in this particular case, by taking h as the parameter A to be varied. Of course we 
could use any of the other parameters appearing in eq. (5.1) but we will limit ourselves 
to this case. At the end of the section we will discuss also eq. (3.1) and then we will also 
use the other parameters appearing in eq. (5.1) as variable parameters A,. We start by 
discussing the model aX h = 0. The potential has a critical line of second-order transitions 
for 6 < and a = 0. At a = b = the model has a tricritical point |p, whereas for 6 > 
there is a line of first-order transitions given by the equation 

a=- (5.2) 

The tricritical point is where the two lines of first- and second-order transitions merge 
together. At any of the first-order transition points one has two phases characterized by 
the following values of (interpreted as the vacuum expectation values, v.e.v.) 

(0)i = 0, (0)2 = ^ (5.3) 
We also define the discontinuity as 

disc (0) = (0)2 - (0)i = ^ (5.4) 

For the following discussion we will assume also a linear dependence of the parameter a 
on the temperature and on the chemical potential around the first-order critical line. We 
could assume an analogous dependence for any of the parameters, but to keep the matter 
simple we will stay with a. We put 

a = aT + jfi — ao (5.5) 

obtaining for the critical line 

aTc + jfic = ao + — (5.6) 
4c 

When we insert the small perturbation h(f) we need (at the lowest order in h) only to know 
how the critical quantities are changed 

aTc + -ffic = ao + \ -rh (5.7) 

4c V 

By noticing that in this case the pressure is given hj p = —V, we get for the specific 
entropy and for the density 

s = -a(02), n = -7(02) (5.8) 
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By taking X = h we can verify eq. (2.10) 



—- disc s H — — disc n 

ah ah 



' ^"^^'^ {-a disc (0^)) + ^ (-7 disc (0^)) 



a\b a dh 

/2c b 



— disc (0) 



(5.9) 



Prom eq. (5.8) and (2.14) we can derive the discontinuity for the energy density 

disc e = {—aTc — 7/Xc) disc (0^) (5.10) 
Then, by differentiating (5.7) 



dh 



1 



a disc {(f)) 



dh 



7 disc {(f)) 



(5.11) 



it is simple matter to verify eq. (2.15). Finally we verify eq. (3.1) by taking h and c as A 
parameters. Using 



disc (0^ 



(5.12) 



we get (in the limit /i — > 0) 



disc 



dTjdc 



dTjdh 



disc {(f)) 



a 4c^ 



-J- 

aV b 



2^ ~ V 2^, 



(5.13) 



All the other relations can be proved in analogous way. 
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